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We have performed ultrasonic measurements on single-crystalline URu 2 Si2 with pulsed mag- 
netic fields, in order to check for possible lattice instabilities due to the hybridized state and 
the hidden-order state of this compound. The elastic constant (Cn-Ci 2 )/2, which is associated 
with a response to the r 3 -type symmetry -breaking (orthorhombic) strain field, shows a three- 
step increase at H > 35 T for H \\ c at low temperatures, where successive meta-magnetic 
transitions are observed in the magnetization. We discovered a new fact that the absolute 
change of the softening of (Cn-Ci 2 )/2 in the temperature dependence is quantitatively recov- 
ered at the suppression of hybridized-electronic state and the hidden order in high-magnetic 
field for H \\ c associated with the successive transitions. The present results suggest that the 
r 3 -type lattice instability, is related to both the emergence of the hybridized electronic state 
and the hidden-order parameter of URu 2 Si 2 . On the other hand, magnetic fields H || [100] 
and [1 10] enhance the softening of (C u -C 'i 2 )/2 in the hidden order phase, while no step-like 
anomaly is observed up to 68.7 T. We discuss the limitation of the localized-electron picture 
for describing these features of URu 2 Si 2 by examination of a crystalline electric field model 
in terms of mean-field theory. 
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'What is the primary order parameter and its ordering vector for the hidden order of 
URu 2 Si 2 ?' is still an open and controversial question and a longstanding issue of condensed 
matter physics. 1 URu 2 Si 2 , which crystallizes in the ThCr 2 Si 2 -type tetragonal structure (space 
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group No. 139, 14/mmm), shows an unknown second-order transition at T a = 17.5 K, known 
as the hidden order (HO) phase, and also a transition into an unconventional superconducting 
state below T c ~ 1.4 K. 1 3 Many theoretical models have been proposed from both local- 
ized and itinerant electron pictures in attempts to identify the order parameter of the HO 
phase. Since 29 Si-NMR and zero-magnetic-field /iSR measurements have detected little (< 
1 Gauss) or no significant internal-magnetic field in the HO phase, 4 6 non-dipolar-type or- 
der parameters, e.g., electric quadrupole (rank 2), hexadecapole (rank 4), r 3 -type magnetic 
octupole (rank 3), or dotriacontapole (rank 5) order parameters, are recently attracting atten- 
tion. 713 Thus far, none of the characteristic signatures of the electric-quadrupole or magnetic - 
octupole order has been identified in resonant X-ray scattering (under magnetic fields) and 
neutron scattering under uniaxial stress. 14 ' 15 On the other hand, a 'nematicity' of the electron 
state in the HO is pointed out, since an in-plane rotational four-fold symmetry breaking is 
observed in the HO by the magnetic-torque measurement by use of cantilever technique. 16 
The interpretation of these experimental results remains controversial. 

URu 2 Si 2 is also considered a heavy-fermion compound since several physical properties 
of this compound exhibit signatures of many-body effects. In particular, the temperature de- 
pendence of the electrical resistivity for / || a and the magnetic susceptibility for H \\ c 
show maxima at around r Amax ~ 70 K and 7^ >max ~ 55 K, respectively, and both quantities 
are strongly suppressed with decreasing temperatures. 117 The whole feature of the electri- 
cal resistivity is similar to the typical Kondo-lattice behavior, whereas the large reduction of 
the c-axis magnetic susceptibility at low temperature seems to be unconventional, suggesting 
the influence of crystalline-electric-field (CEF) effects and/or large antiferromagnetic corre- 
lations. 7111218 However, there is no consensus of the mechanism behind the formation of 
the low-temperature heavy-electron state of this system, even on the CEF 5f level schemes 
assumed in the local limit. 

On the other hand, when applying a magnetic field over 35 T along c-axis, three- 
successive transitions have been observed in resistivity, specific-heat, magnetization, and 
longitudinal-ultrasonic-velocity measurements. 17 ' 19-22 These cascade transitions occur in the 
vicinity of the upper phase boundary of the HO phase and are associated with meta-magnetic- 
like increases in the c-axis magnetization. 19 23 Intriguingly, the crossover temperatures of 
T p>max and 7^ >max are gradually reduced with increasing magnetic field and vanishe in the 
vicinity of the HO phase boundary in the magnetic-field-temperature phase diagram. 17 This 
fact suggests that a drastic change from the heavily hybridized to little or non-hybridized 
electronic state is expected with the collapse of the HO. Here, one of the outstanding issues 
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is whether any signature of symmetry breaking can be observed within the cascade-transition 
region in order to clarify the characteristics of the HO and multiple phases in high magnetic 
field. 

Ultrasonic measurements can sensitively detect symmetry-breaking lattice instabilities 
and electric quadrupole responses in the single crystals. Early studies employing ultrasonic 
measurements on URu 2 Si 2 uncovered a characteristic decrease (softening) of the elastic con- 
stant Cn and (Cn-Ci 2 )/2 with decreasing temperature. 24 ' 25 The softening of C n had been 
explained by Griineisen coupling due to many-body effects, which appears in the bulk mod- 
ulus of typical heavy-fermion compounds. This phenomena is the so-called Kondo-volume 
collapse. 26 On the other hand, a study of elastic properties of URu 2 Si 2 in comparison with 
the non-5/ contribution of ThRu 2 Si 2 implies that the softening of (Cn-Ci 2 )/2 in URu 2 Si 2 
originates from uranium's 5/ electrons. 27 The (C n -Ci 2 )/2 mode, which corresponds to or- 
thorhombic strain that can conserve volume and breaks T 3 symmetry, is the only transverse 
mode which exhibits softening in this compound. Such a mode- selective elastic response re- 
minds us of a T 3 -type crystal deformation due to orbital fluctuations of c-f hybridized bands 
(band Jahn-Teller effect) as an origin of softening rather than a CEF effect with a pseudo de- 
generate ground state, which potentially causes softening in other transverse modes. They are, 
however, not fully understood by both itinerant and localized 5/-electrons pictures, since the 
c-f hybridized band structure and possible contribution of the band for the elastic responses 
have not been clarified, thus far. 

The present study addresses the above issues by means of ultrasonic measurements 
on single-crystalline URu 2 Si 2 with pulsed-magnetic fields up to 68.7 T at the Hochfeld- 
Magnetlabor Dresden to check for the elastic response at high magnetic fields. In the present 
paper, we report on the magnetic-field dependence of the trans verse- ultrasonic mode (Cn- 
Ci 2 )/2 for H || c and also for H ± c. This provides new information about the relationship 
between the HO and symmetry-breaking lattice instability. So far, there are only a few re- 
ports concerning ultrasonic studies of the cascade transitions in the high field region using 
the longitudinal Cn and C 33 modes, which correspond to strain that changes volume while 
preserving the tetragonal symmetry. 20 - 22 

The sample used was a single crystal of URu 2 Si 2 with the dimensions of 3.8 x 1.8 x 
1.2 mm 3 ([1 10]-[1 10]-[001]). LiNb0 3 transducers with thickness of 40 yum are fixed on the 
mirror-polished [1 10] surfaces and used for generation and detection of ultrasound. The sound 
velocity was measured by use of the phase-comparator method by recording the in-phase and 
quadrature signals with repetition rate of 20 kHz at fixed frequencies 255.6 MHz and 155.3 
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MHz for// || c and H ± c, respectively. In the present measurements, an absolute value of the 
sound velocity for propagation A: || [110] and polarization « || [110] wasv = 2511 m/s at 4.2 K, 
which results in an absolute value of the elastic constant (C n -Ci 2 )/2 = 6.453xl0 10 J/m 3 with 
the density p = 10.01 g/cm 3 of URu 2 Si 2 . The pulsed-magnetic fields had a duration of ~25 
ms with maximum field of 46.8 T for H || c and a duration of ~150 ms with maximum field 
of 68.7 T for H ± c. As shown below, the signal-to-noise ratio of the present measurement 
is small for all measurements, indicating that the measurement conditions were maintained 
during the sample rotations. In particular, the sample was manually rotated together with a 
non-magnetic holder without removing the transducers and signal lines from the sample for 
the in-plane magnetic field while rotating from H || [100] to [1 10]. 

Figure 1 shows the magnetic-field dependence of (Cn-Ci 2 )/2 and of the ultrasonic- 
attenuation coefficient Aa at 1.5 K. (Cn-Ci 2 )/2 incleases drastically with H \\ c, but de- 
creases moderately for H ± c. This anisotropy shows that the origin of the elastic softening 
observed in the HO state is being suppressed for H \\ c while being enhanced for H ± c. 
For H || c, we can see three anomalies in the successive-transition region, 35 T < H < 39 
T, as indicated by the downward pointing arrows in upper panel of Fig. 1 and inset (a). At 
4.2 K, only one kink, which corresponds to the upper HO-phase boundary, can be identified. 
These transition fields are consistent with previous reports of pulsed- and static-magnetic 
field measurements. 17 ' 23 ' 28 This indicates that an isothermal condition was maintained during 
the pulse in the present measurement, i.e., a magneto-caloric effect (magnetic refrigeration) 
in the vicinity of the phase boundary is not occurring, while such effects easily occur when 
the pulse duration is too short. Indeed, the elastic anomalies shift to lower fields and broaden 
at 4.2 K as shown in inset (a) of Fig. 1 . On the other hand, no elastic anomaly is observed for 
H J. c up to 68.7 T as shown in inset (b) of Fig. 1. The absolute value of the elastic constant 
after applying pulsed magnetic fields (at H = 0) always returns to the value observed before 
the pulse. This shows that there is no eddy-current heating even for the large cross section 
(~6.8 mm 2 ) and short pulse-duration for the magnetic field aligned along c. Incidentally, no 
clear elastic anomaly was detected in the HO phase for H || c at around H* ~ 22 T, where a 
reconstruction of the Fermi surface is observed from electrical transport measurements using 
high quality single crystal. 29 30 A small hysteresis-like loop of uncertain origin is, however, 
observed in the intermediate magnetic-field region. 

The in-plane magnetic field dependence for H || [100] and H || [110] in the HO phase at 
1.5 and 3.9 K is shown in inset (b) of Fig. 1. The slight difference of 4xlO -5 in the response of 
(Cn-Ci 2 )/2 can be recognized between the data for H || [100] and [1 10] in the high magnetic- 
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Fig. 1. (Color online) Magnetic-field dependence of elastic constant (C\\-C\i)l2 (upper panel) and sound- 
attenuation change Aa (lower panel) for H \\ c and H X c at 1.5 K. Inset (a) shows elastic anomalies in the 
vicinity of the cascade transitions at 1.5 and 4.2 K. Inset (b) shows (C\\-C\2)j2 vs. magnetic field for H || [100] 
and H || [110] at 1.5 and 3.9 K The data, except the inset (a), show both up- and down-sweep of the magnetic 
fields. 



field region. The difference is, however, within the measurement accuracy of ~ 5 x 10~ 5 from 
the noise signal. Thus, we conclude that there is no evidence of in-plane anisotropy. Later, we 
indeed will show that the expected effect is too small to be observable in our measurement. 

In Fig. 2, we present the temperature and magnetic-field dependence of the elastic re- 
sponse for two individual single crystals with similar residual resistivity ratio (RRR = 
P(i.5K)/P(3oqk))- Sample#l has RRR ~ 10 and sample#2 has RRR ~ 19. These two data sets 
are plotted with the same ordinate; the upper axis corresponding to the blue data shows the 
temperature dependence at T, the lower axis corresponding to the data shows the magnetic- 
field dependence at 1.5 K. Here, we notice that the elastic-constant change (softening of 
~0.7%) in the temperature dependence below the local maximum at ~ 120 K is comparable 
to the increase up to H ~ 45 T in the field dependence. Here, the maximum in the elastic 
constant eventually occurs as a balance between the phonon-background and the softening 
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Fig. 2. (Color online) Comparison of elastic constant change in the temperature dependence (blue plots) and 
magnetic field dependence for H \\c (red plots). 



of strain-susceptibility due to 5 /-electron effects in the temperature dependence, 27 while the 
c-axis magnetization shows a maximum at around T Xtmsa ~ 55 K. It should be mentioned that 
the present temperature dependence is different from the previously reported data, 25 which 
shows the maximum at ~ 70 K and also the softening ~ 0.2% for (C\\-Cu)l2. We consider 
these differences are due to the recent improvements of the measurement conditions, e.g., 
using mirror-polished sample, appropriate adhesive agents, and higher frequency (higher di- 
rectional quality) of the ultrasonic wave. It should be noted that the c-axis magnetization 
also responds similarly to temperature and magnetic field changes. So we point out that 
there are considerable similarities between the reduction of c-axis magnetization and elas- 
tic constant {C\\-Cyi)l2 in the temperature dependence, and also the different tendencies of 
the elastic constant between H \\ c and H _L c in the magnetic field dependence and the 
strong anisotropy in the c- and a-axis magnetization in URu 2 Si2- These similarities strongly 
imply that these phenomena will share a common root cause. The origin must be related to 
the low-temperature electronic state of this compound and also an anisotropic nature of the 
HO parameter, e.g., the higher-rank multipole order or the c-f hybridized-band instability. 
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Fig. 3. (Color online) (a): Calculated (uniform) magnetic- and quadrupole-susceptibilities vs. T, and (b): 
calculated (uniform) magnetic- and quadrupole-susceptibilities vs. H for H || [001], [110] and [100] at 1.5 K up 
to 70 T using the MF theory with a CEF scheme as described in the text. 



As described above, since the hybridization effect will be dominant in the low- 
temperature and low-magnetic field region of this compound, it is unreasonable to repro- 
duce the bulk susceptibilities only based on the localized electron model. Nevertheless, some 
features can be explained by assuming an appropriate CEF level scheme which mimics the 
hybridization effect, including the band Jahn-Teller effect, even in the localized 5/-electron 
picture. 7 ' 11 In this case, it can be expected that the recovery of the elastic constant (Cn-Ci2)/2 
around the cascade transitions in the high magnetic field region is accompanied by a drastic 
change of the CEF ground state. In the present paper, we dare to introduce a CEF model 
in terms of mean-field (MF) theory assuming an antiferro-hexadecapole (AFH) order which 
reproduce the temperature and magnetic field dependence of (C\\-Cu)l2 in both the normal 
and hybridized states, in order to display the limitations of the localized-electron picture for 
describing these responses of (Cn-C l2 )/2. 

Figure 3 (a) shows CEF calculations of uniform susceptibilities vs. T, and Fig. 3 (b) shows 
the magnetic-field dependence of these uniform susceptibilities at 1 .5 K. These uniform mag- 
netic and quadrupole susceptibilities are obtained by numerical differentiation with the tiny 
magnetic and strain field, respectively. The present CEF level scheme is; F] (0 K)-r 2 (60 
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K)T 3 (178 K)-Y\ l> (491 K), for the Hamiltonian, // CEF = - Z m ,„ B n m O n m (I, m = 0,2, 4, 6), with 
B° 2 = —12 K, fl° = —0.43 K, B\ = —3.2 K, fl° = —0.011 K, and B\ = 0.053 K, and the 
MF parameters (as will be described later). Note that the present CEF model 31 has a similar 
level scheme as originally proposed for PrRu 2 Si 2 . 32 Based on the localized-electron picture, 
the temperature and magnetic-field dependence of (Cn-Ci2)/2 is understood as a quadrupole 
susceptibility ~xr3 with the r 3 -symmetry quadrupole operator O v = -^0\, which is defined 
as O v = (7 2 + J 2 _)/2. The present CEF calculation can reproduce the following experimental 
facts; an anisotropy between c- and a-axis magnetizations at high temperature region T > 100 
K, a softening of {C n -C i 2 )/2 in the temperature dependence, and also an anisotropic response 
of (Cn-Ci 2 )/2 and magnetization for H \\ c and H ± c. Needless to say, other CEF models, 
which have been proposed thus far, can not reproduce the anisotropic elastic response, e.g., 
Santini's 5f 2 model, 7 r 4 -r ( 1 1) (44 K)-r 2 (111 K), and Garatanu's 5f 2 model, 18 r^-I^ (404 
K)-r 2 (1076 K), give rise to softening for C 44 or Nieuwenhuys's model, 11 r ( J 1) -r 2 (50 K)-r (2) 
(171 K), does not even reproduce the softening of (Cn-Ci 2 )/2. 27 

Regarding the possibilities of the AFH order scenario on the localized-electron picture, 
which has recently been proposed by Kusunose and Harima, 11 the present CEF scheme, which 
includes a r 2g -symmetry xy(x 2 - j 2 )-type hexadecapole H" as an active multipole, does not 
forbid the AFH order. By using the mean-field (MF) theory of the //"-type AFH order with 
the CEF model, the effective Hamiltonian is written as, H^ F = ± [D h H z (i)H"(j) + 
D d J z (i)J z (j)] + Xii^cEF - 8jI^bH zZiJ z (i). Here, D h and D d are coupling constants for the 
inter-site electric hexadecapole and magnetic dipole interactions, respectively. The additional 
depression of (Cn-Ci 2 )/2 just below T Q in the temperature dependence can also be repro- 
duced by the quadrupole susceptibility ~Xr3 assuming the MF parameters as, Dd = 2.093 
K, Dh = 0.0433 K. On the other hand, the present MF model can not reproduce a tempera- 
ture dependence of the elastic constant C^, where the calculated quadrupole susceptibility of 
Oxy = (Jl - J-)/2i decreases with decreasing temperature below T D while the experimental 
result of C66 exhibits an opposite temperature dependence in the HO phase. 24 ' 25 Moreover, 
the magnetic-field dependence of ~xn for H \\ [001] does not show any anomaly at the sup- 
pression of the HO, while the calculated c-axis magnetic susceptibility shows a jump at H\. 
These discrepancies simply show the limitations of the present localized-electron model to 
describe the elastic responses. 

After all, it is undeniable whether the softening of C\\ is caused by the change of 
bulk modulus Cb due to Kondo effect or a softening of (Cn-Ci 2 )/2, since the C\\ mode in 
the tetragonal symmetry includes both contributions, i.e., it is decomposed as C\\ = 3C B - 
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C„+(Cii-Ci 2 )/2-2Ci3 with bulk modulus C B = (2C n + 2Ci 2 + AC n + C 33 )/9 and tetragonal 
strain mode C u = (C n + C n - 4Ci 3 + 2C 33 )/6. In order to estimate the contribution of C B 
on C n , further ultrasonic measurements will be needed to check the longitudinal mode, with 
k || u |] [110], C L no = 3Cb-C„+C 66 -2Ci 3 , and the other transverse modes C 44 and C 66 for 
comparison. 

Finally, we will focus on the in-plane magnetic-field dependence at 1 .5 K, where no phase 
transition is found up to 68.7 T (Fig. 3(b)). The difference of -^ r3 between H || [100] and 
[1 10] in the calculated results in the high magnetic field region is mainly attributable to the 
CEF effect with or without MF interactions. Despite this, it can be expected that the bulk 
susceptibility must be affected by induced lower-rank multipoles when applying high mag- 
netic fields H ± c in the AFH ordered phase. For example, a magnetic dipole and an electric 
quadrupole should be induced via Ginzburg-Landau coupling. 31 When the change of cal- 
culated susceptibility is scaled to the actual elastic constant change from T to 60 T, the 
contribution of the difference between H || [100] and [110] is estimated as ~ 0.5 x 10~ 5 at 60 
T. The magnitude of the modulation due to MF interactions accounts for less than 1 % of this 
small anisotropy. As described above, even the possible difference between H || [100] and 
[110] due to the CEF effect could not be distinguished within the present margin of the error 
~ 5x 10~ 5 . 

In conclusion, we have performed, for the first time, the ultrasonic measurement of (C n - 
Cn)/2 in the HO phase of URu 2 Si 2 for H \\ c and H ± c under pulsed magnetic fields up to 
68.7 T. Elastic anomalies of (Cn-Ci 2 )/2, which is a symmetry-breaking (volume-conserving) 
ultrasonic mode, is observed at around the cascade transition region 35 T < H < 39 T. For 
H ± c, we found that there is no clear difference between H || [100] and [110] within the 
present measurement accuracy of ~ 5 x 10~ 5 . The present experimental results suggest that 
the suppression of the HO in high magnetic field is accompanied by not only the recovery 
from the heavy-fermion state, which causes a reduction of the c-axis magnetic moment at 
low temperature, but also the reduction of the T 3 lattice instability. We conclude that the low- 
temperature electronic state of URu 2 Si 2 has a T 3 -type (jc 2 -y 2 -type) lattice instability which is 
related to the origin of the hybridized electronic state and the HO, such as higher-multipole 
fluctuations and/or band Jahn-Teller effects with aT 3 symmetry. Thus, an understanding of the 
55-120 K energy-scale phenomena due to strong hybridization based on the itinerant electron 
pictures, as seen in the magnetization and elastic constant (Cn-Ci 2 )/2, should be refocused 
to elucidating the true nature of the HO. 
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